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Abstract 

We study the problem of finding the worst-case joint distribution of a 
set of risk factors given prescribed multivariate marginals and a nonlin¬ 
ear loss function. We show that when the risk measure is CVaR, and the 
distributions are discretized, the problem can be conveniently solved using 
linear programming technique. The method has applications to any situa¬ 
tion where marginals are provided, and bounds need to be determined on 
total portfolio risk. This arises in many financial contexts, including pric¬ 
ing and risk management of exotic options, analysis of structured finance 
instruments, and aggregation of portfolio risk across risk types. Applica¬ 
tions to counterparty credit risk are emphasized, and they include assessing 
wrong-way risk in the credit valuation adjustment, and counterparty credit 
risk measurement. Lastly a detailed application of the algorithm for coun¬ 
terparty risk measurement to a real portfolio case is also presented in this 
paper. 
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1 Introduction 


In recent years counterparty credit risk management has become an increasingly 
important topic for both regulators and participants in over-the-counter deriva¬ 
tives markets. Even before the hnancial crisis, the Counterparty Risk Management 
Policy Group noted that counterparty risk is “probably the single most impor¬ 
tant variable in determining whether and with what speed hnancial disturbances 
become hnancial shocks, with potential systemic traits” (CRMPG [2005]). This 


concern over counterparty credit risk as a source of systemic stress has been re- 


hected in the historical developments in the Basel Capital Accords (BCBS 
BCBS [ 2011 ] see also Section 3 below). 


2006 


Counterparty Credit Risk (CCR) is dehned as the risk of loss due to default or the 
change in creditworthiness of a counterparty before the hnal settlement of the cash 
hows of a contract. An examination of the problems of measuring and managing 
this risk reveals a number of key features. First, the risk is bilateral, and current 
exposure can he either with the institution or its counterparty. Second, the evalu¬ 
ation of exposure must be done at the portfolio level, and must take into account 
relevant credit mitigation arrangements, such as netting and the posting of col¬ 
lateral, which may be in place. Third, the exposure is stochastic in nature, and 
contingent on current market risk factors, as well as the creditworthiness of the 
counterparty, and credit mitigation. In addition, possible dependence between 
credit risk and exposure, known as wrong-way risk, is also an important mod¬ 
elling consideration]^ Finally, the problem is computationally highly intractable. 
To calculate risk measures for counterparty credit risk, one requires the joint 
distribution of all market risk factors ahecting the portfolio of (possibly tens of 
thousands of) contracts with the counterparty, as well as the creditworthiness of 
both counterparties, and values of collateral instruments posted. It is nearly im¬ 
possible to estimate this joint distribution accurately. Consequently, we are faced 
with a problem of risk management under uncertainty, where at least part of the 
probability distribution needed to evaluate the risk measure is unknown. 


Fortunately, we do have partial information to aid in the calculation of coun¬ 
terparty credit risk. Most hnancial institutions have in place models for sim¬ 
ulating the joint distribution of counterparty exposures, created (for example) 
for the purpose of enforcing exposure limits. Additionally, internal models for 
assessing default probabilities, and credit models (both internal and regulatory) 
for assessing the joint distribution of counterparty defaults are also available at 

^Since in general the risk is bilateral, in the case of pricing contracts subject to counterparty 
credit risk (i.e. computing the credit valuation adjustment), the creditworthiness of both parties 
to the contract is relevant. In this paper, we take a unilateral perspective, focusing exclusively 
on the creditworthiness of the counterparty. 


2 














our disposal. We can view the situation as one where we are given the (multi¬ 
dimensional) marginal distributions of certain risk factors, and need to evaluate 
portfolio risk for a loss variable that depends on their joint distribution. The 
Basel Accord (BCBS [2006]) has employed a simple adjustment based on the “al¬ 
pha multiplier” to address this problem. A stress-testing approach, employing 
different copulas and financially relevant “directions” for dependence between the 


market and credit factors is presented in Garcia-Cespedes et ah 2010 and Rosen 
and Saunders [2010 . This method allows for a computationally efficient evalua¬ 


tion of counterparty credit risk, as it leverages pre-computed portfolio exposure 

simulations El 

In this paper we investigate the problem of determining the worst-case joint dis¬ 
tribution, i.e. the distribution that has the given marginals, and produces the 
highest risk measure E] This approach is motivated by a desire to have conserva¬ 
tive measures of risk, as well as to provide a standard of comparison against which 
other methods be evaluated. While in this paper we focus on the application to 
counterparty credit risk, as mentioned earlier, the problem formulation is com¬ 
pletely general, and can be applied to other situations in which marginals for risk 
factors are known, but the joint distribution is unknown. Finally, we note that we 
work with Conditional Value-at-Risk (CVaR), rather than Value-at-Risk (VaR), 
which is the risk measure that currently determines regulatory capital charges 


for counterparty credit risk in the Basel Accords (BCBS [2006], BCBS [2011]). 


The motivation for this choice is twofold. First, it yields a computationally more 
tractable optimization problem for the worst-case joint distribution, which can be 
solved using a linear programming technique. Although not directly in the con¬ 
text of CCR, the Basel Committee is also considering the replacement of VaR with 
CVaR as the risk measure for determining capital requirements for the trading 


book (BCBS [2012]). 


Model uncertainty, and problems with given marginal distributions or partial in¬ 
formation have been studied in many financial contexts. One example is the 
pricing of exotic options, where no-arbitrage bounds may be derived based on 


observed prices of liquid instruments. Related studies include 

Bertsimas and 

Popesci 

1 2002 , 

Hobson et al. 

2005a , 

Hobson et al 

2005b , 

Laurence and Wang 

2004 , 

jaurence and Wang |2005|, and 

Chen et al. 

2008 , for exotic options writ- 


ten on multiple assets {Si, ..., St) observed at the same time T. The approach 


closest to the one we take in this paper is that of Beiglbock et al. 2011 


m 


^Generally speaking, the computational cost of the algorithms is dominated by the time 
required for evaluating portfolio exposures - which involves pricing thousands of derivative con¬ 
tracts under at least a few thousand scenarios at multiple time points - rather than from the 
simulation of portfolio credit risk models. 

^Since the marginals are specified, the problem is equivalent to finding the worst-case joint 
distribution with the prescribed (possibly multi-variate) marginals. 
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which the marginals ('h(S'f),..., 'h(S'J)) are assumed to be given, and an inhnite- 
dimensional linear programming technique is employed to derive price bounds. 
There is also a large literature on deriving bounds for joint distributions with given 
marginals, and corresponding VaR bounds. For a recent survey, see|Puccetti and 


Riischendorf 2012 . The problems considered in this paper can be characterized 
by two important aspects; (i) we use an alternative risk measure (CVaR), which 
are provided with multivariate (non-overlapping) marginal distributions, and (ii) 
we have losses that are a non-linear (and non-standard) function of the underlying 
risk factors. Glasserman and Xu |2014| present a similar empirical approach to 
estimating worst-case error in a range of risk management problems; in addition 
to discussing the theoretical aspects of this problem, one of the primary goals 
of our paper is addressing the numerical challenges arising from worst-case joint 
distribution problem. 


Haase et al. 2010 propose a model-free method for a bilateral credit valuation 


adjustment; in particular their proposed approach does not rely on any specific 
model for the joint evolution of the underlying risk factors. [Talay and Zhang 


2002 treat model risk as a stochastic differential game between the trader and the 


market, and prove that the corresponding value function is the viscosity solution of 


the corresponding Isaacs equation. Avellaneda et al. 1995 , Denis et al. 2011 and 


Denis and Martini 2006 consider pricing under model uncertainty in a diffusion 


context. Recent works on risk measures under model uncertainty include Kervarec 


2008 and Bion-Nadal and Kervarec 2012 


The remainder of the paper is structured as follows. Section 2 frames the problem 
of hnding worst-case joint distributions for risk factors with given marginals, and 
show how this can be reduced to a linear programming problem when the risk 
measure is given by CVaR and the distributions are treated as discrete. Section 
3 outlines the application of this general approach to counterparty credit risk in 
the context of the model underlying the CCR capital charge in the Basel Accord. 
In section 4 a nontrivially numerical example using a real portfolio is provided, 
and section 5 presents conclusions and directions for future research. 


2 Worst-Case Joint Distribution Problem 


Let Y and Z be two vectors of risk factors. We assume that the multi-dimensional 
marginals of Y and Z, denoted by Fyi^) and Fzi^) respectively, are known, 
but that the joint distribution of (V, Z) is unknown (Note: in the context of 
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counterparty credit risk management discussed in the next section, Y and Z will 
be vectors of market and systematic credit factors respectively). Portfolio losses 
are defined to be L = L{Y, Z), where in general this function can be non-linear. 
We are interested in determining the joint distribution of {Y, Z) that maximizes 
a given risk measure p: 

max piLiY Z)) (1) 

mrYz) 

where diFy, Fz) is the Frechet class of all possible joint distributions of {Y, Z) 
matching the previously dehned marginal distributions Fy and Fz- More ex¬ 
plicitly, for any joint distribution Fyz € '^{Fy,Fz) we have Ily{Fyz} = Fy and 
^z{Fyz} = Fz, where If {.} denote the projections that take the joint distribution 
to its (multi-variate) marginals. While we are mainly interested in the application 
aspect of this, we note that bounds on instrument prices can be derived within 
the above formulation by taking the risk measure to be the expectation operator. 

It is well known that given a time horizon and conhdence level a, Value-at-Risk 
(VaR) is dehned as the a-percentile of the loss distribution over the specihed time 
horizon. The shortcomings of VaR as a risk measure have been much discussed in 
the literature. An alternative measure that addresses many of these shortcomings 
is Conditional Value at Risk (CVaR), also known as tail VaR or Expected Shortfall. 
If the loss distribution is continuous, CVaR is the expected loss given that losses 
exceed VaR. More formally, we have the following dehnition of CVaR. 

Definition 2.1. For the confidence level a and loss random variable L, the Con¬ 
ditional Value at Risk at level a is defined by 

1 C 

CVaR„(L) = / VaR^{L)df 

^ rv 


2008 


We will use the following result, which is part of a theorem from Schied 
(translated into our notation). Here L is regarded as a random variable dehned 
on a probability space (H, B, F). 


Theorem 2.1. CVaRQ(L) can be represented as 


CVaRo = sup E(e[L] 

GePc 

where Qa is the set of all probability measures G ]F|^ whose density dG/dF is 
F-a.s. bounded by 1 /{1 — a). 

<C F means G is absolutely continuous with respect to F, i.e. for any b G B that F(B) = 0 
we have G{B) = 0. 
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Applying the above result, with p = CVaRa, the worst-case joint distribution 
problem stated in Q can be conveniently reformulated as: 

sup Eg[L] (2) 

F,Ged{FY,Fz) 

U,{F} = Fy 

n,{F} = Fz 

dG 1 

—— ^ - a.s. 

dF 1 — a 

Note that the hnal constraint assumes explicitly that the corresponding density 
exists. 

In many practical cases the marginal distributions will be discrete, either due 
to a modelling choice, or because they arise from the simulation of separate 
continuous models for Y and Z. In this case, the marginal distributions can 
be represented by FyiY = ym) = Pm,rn = I,M, and Fz{Z = Zn) = Qn, 
n = 1,... ,N. Any joint distribution of {Y, Z) is then specihed by the quantities 
FyziY = Hm-iZ = Zn) = V'mn, and the worst-case CVaR optimization problem 
above can be further simplihed to: 


max ■ 

1 — Of 

n,m 

pmn 

(3) 

^ ^ Pm ^ 1; • 

..,M 


^ ^ Qn ^ I? • 

..,N 


^ ^ P^mn 1 ^ 




n,m 

0 ^ Pmn ^ '4^mm 


Evidently this is a well-dehned linear programming problem, and has the general 
form of a mass transportation problem with multiple constraints. Note that since 
the sum of each marginal distribution is equal to one, we do not have to include the 
additional constraint that the total mass of p equals one. Excluding the bounds, 
this LP has 2mn variables and constraints. Consequently, the above 

formulation can lead to very large linear programs. In the numerical examples pre¬ 
sented later in this paper, we employ marginal distributions with market scenarios 
and credit scenarios ranging from 1,000 to 5,000, yielding joint distributions of 
0 ( 10 ^). 

Specialized algorithms for linear programs that take advantage of the structure of 
the transportation problem may well be required for problems dehned by marginal 
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distributions with larger numbers of scenarios. The reader if referred to Rachev 


and Riischendorf 1998 for more information on various types of Transportation 


Problem. 


3 Wrong-way Risk and Counterparty Credit Risk 


The internal ratings based approach in the Basel Accord ( |BCBS 2006 ) provides a 
formula for the charge for counterparty credit risk capital for a given counterparty 
that is based on four numerical inputs: the probability of default (PD), exposure 
at default (EAD), loss given default (LGD) and maturity (M). 


Capital = EAD ■ LGD 


$ 


^-\PD) + ^ - $- 1 ( 0 . 999 ) 
V^-p 


MA(M,PD) (4) 


Here $ is the cumulative distribution function of a standard normal random vari¬ 
able, and MA is a maturity adjustment (see BCBS 2006 The probability of 


default is estimated based on an internal rating system, while the LGD is the 
estimate of a downturn LGD for the counterparty based on an internal model. 
Another parameter appearing in the formula, the correlation (p), is essentially 
determined as a function of the probability of default. 

The exposure at default in the above formula is a constant. However, as noted 
above, counterparty exposures are inherently stochastic in nature, and potentially 
correlated with counterparty defaults (thus giving rise to wrong-way risk). The 
Basel accord circumvents this issue by setting EAD = a x Effective EPE, where 
Effective EPE is a functional of a given simulation of potential future exposures 
(see BCBS 2006 , De Prisco and Rosen 2005 or Garcia-Cespedes et al. 2010 


for detailed discussions). The multiplier a defaults to a value of 1.4; however 
it can be reduced through the use of internal models (subject to a floor of 1.2). 
Using internal models, a portfolio’s alpha is dehned as the ratio of CCR economic 
capital from a joint simulation of market and credit risk factors (E'G^°*“^) and 
the economic capital when counterparty exposures are deterministic and equal to 
expected positive exposure]^ 

jp/^Total 

(5) 


a = 




^In the most recent version of the charge, exposure at default may be reduced by current 
CVA, and the maturity adjustment may be omitted, if migration is accounted for in the CVA 
capital charge. See BCBS [2006] for details. 

^Expected positive exposure (EPE) is the average of potential future exposure, where av¬ 


eraging is done over time and across all exposure scenarios. See BCBS [2006][ De Prisco and 
Rosen 2005 or Garcia-Cespedes et al. 2010 for detailed expressions of this. 
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The numerator of a is economic capital based on a full joint simulation of all 
market and credit risk factors (i.e. exposures are treated as being stochastic, and 
they are not treated as independent of the credit factors). The denominator is 
economic capital calculated using the Basel credit model with all counterparty 
exposures treated as constant and equal to EPE. For infinitely granular portfolios 
in which PFEs are independent of each other and of default events, one can assume 
that exposures are deterministic and given by the EPE. Calculating a tells us how 
far we are from such an ideal case. 


3.1 Worst-Case Joint Distribution in the Basel Credit Model 


In this section we demonstrate how the worst-case joint distribution problem can 
be appliedto the Basel portfolio credit risk model for the purpose of calculating 
the worst-case alpha multiplier. 

In order to calculate the total portfolio loss, we have to determine whether each 
of the counterparties in the portfolio has defaulted or not. To do so, we define the 
creditworthiness index of each counterparty k,l ^ k ^ K, using a single factor 
Gaussian copula a^ 

CWIfc = y/pk ■ Z yjl — Pk ■ €k (6) 

where Z and are independent standard normal random variables and pk is the 
factor loading giving the sensitivity of counterparty k to the systematic factor Z. 
If PDfc is the default probability of counterparty fc, then that counterparty will 
default if: 

CWIfc < ^-\FDk) 


Assuming that we have M < oo market scenarios in total, if ykm is the exposure 
to counterparty k under market scenario m, the total loss under each market 
scenario is: 

K 

L„, = Y,ykm-l {CWlk < $-'(PDfc)} (7) 

k=l 


Below we focus on the co-dependence between the market factors Y and the 
systematic credit factor Z. In particular, we assume that the market factors 
Y and the idiosyncratic credit risk factors Sk are independent. This amounts 
to assuming that there is systematic wrong-way risk, but no idiosyncratic wrong¬ 


way risk (see Garcia-Cespedes et ah 2010 for a discussion). Define the systematic 


^In principle we can introduce a fat-tailed copula in place of the Gaussian copula, but with 
a considerable increase in the resulting computational requirement. 







losses under market scenario m to be: 


a/1 - Pk 


( 8 ) 


K 


Lm{Z) / ^ Vh 




k=l 


with probability P(y = ym) = Pm- Next we discretize the systematic credit factor 
Z using N points and define Lmn as: 


K 

Lmni^Z^ ^ ^ ykm^ 
k=l 


‘1’ (PDfc) — y/f^ ■ Z„ 

a/1 - Pk 


F{Z = Zn) = Qn for 1 < n < 

where L^n represents the losses under market scenario m, 1 < m < M, and credit 
scenario n, 1 < n < A^. 


In hnding the worst-case joint distribution, we focus on systematic losses, and sys¬ 
tematic wrong-way risk, and consequently we need only discretize the systematic 
credit factor Z. We employ a naive discretization of its standard normal marginal: 


Fz{Z = Zn) = qn = ^{Zn) - ^{Zn-l) j = 1, . . . , N 


where Zq = —oo and z^^i = oo. In the implementation stage in this paper, we 
set N = 1000, and take Zj to be eqnally spaced points in the interval [—5, 5]. This 
enables us to consider the entire portfolio loss distribntion nnder the worst-case 
joint distribution. In calculating risk at a particular conhdence level, there is 
potentially still much scope for improvement over our strategy by choosing a finer 
discretization of Z in the left tai0 


For a given confidence level a, the worst-case joint distribntion of market and 
credit factors, 'ipmn, = 1, • • •, M, n = 1,... ,N can be obtained by solving the 
LP stated in ([^. Having fonnd the discretized worst-case joint distribution, we 
can simulate from the full (not just systematic) credit loss distribution using the 
following algorithm in order to generate portfolio losses: 


1. Simulate a random market scenario m and credit state N from the discrete 
worst-case joint distribution 

2. Simnlate the creditworthiness index of each connterparty. Snpposing that 
Zn is the credit state for the systematic credit factor from Step 1, simulate Z 
from the distribution of a standard normal random variable conditioned to 

®Although we have used an evenly spaced grid for discretizing Z, importance sampling tech¬ 
niques can be utilized to better capture the behaviour of worst-case joint distribution in the left 
tail. 


9 








be in (zn-i,Zn)- Then generate K i.i.d. standard normal random variables 
Ek, and determine the creditworthiness indicators for each counterparty using 
equation ([^. 

3. Calculate the portfolio loss for the current market/credit scenario: using the 
above simulated creditworthiness indices and the given default probabilities 
and asset correlations, calculate either systematic credit losses using ([^ or 
total credit losses using ([^. 

4 Application to Counterparty Credit Risk 


In this section we consider the use of the worst-case joint distribution problem 
to calculate an upper bound on the alpha multiplier for counterparty credit risk 
using a real-world portfolio of a large financial institution. The portfolio con¬ 
sists of over-the-counter derivatives with a wide range of counterparties, and is 
highly sensitive to many risk factors, including interest rates and exchange rates. 
Results calculated using the worst-case joint distribution are compared to those 
using the stress-testing algorithm correlating the systematic credit factor to to¬ 
tal portfolio exposure, as described in Garcia-Cespedes et ah 2010| and |R 


tosen 


and Saunders 2010 . More specihcally, we begin by solving the worst-case CVaR 


linear program (|^ for a given, pre-computed set of exposure scenarios]^ and the 
discretization of the (systematic) credit factor in the single factor Gaussian cop¬ 
ula credit model described above. We then simulate the full model based on the 
resulting joint distribution, under the assumption of no idiosyncratic wrong-way 
risk (so that the market factors and the idiosyncratic credit risk factors remain 
independent of each other). 


The market scenarios are derived from a standard Monte-Garlo simulation of 
portfolio exposures, so that we have: 


Pv(R = 2/m) = Pm = ^ i = l,...,M (9) 

In the coming section we will look at top counterparties with respect to total 
portfolio exposure and some of their exposure characteristics. 

^Exposures are single-step EPEs based on a multi-step simulation using a model that assumes 
mean reversion for the underlying stochastic factors. 
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4.1 Portfolio Characteristics 


The analysis presented in this section is based on a large portfolio of over-the- 
connter derivatives inclnding positions in interest rate swaps and credit defanlt 
swaps with approximately 4,800 connterparties. We focns on two cases, the largest 
220 and largest 410 connterparties as ranked by exposnre (EPE); these two cases 
acconnt for more than 95% and 99% of total portfolio exposnre respectively. 

Figures [T] and [^present exposure concentration reports, giving the number of ef¬ 
fective counterparties among the largest 220 and 410 counterparties respectively}^ 
The effective number of counterparties for the entire portfolio in shown in Figure 
3. As can be seen in these hgures the choice of largest 220 and 410 counterparties 
is justihed as the number of effective counterparties for the entire portfolio is 31 
in each case. 


Effective Counterparty Number 



Figure 1: Effective number of counterparties for the largest 220 counterparties. 

The exposure simulation uses M = 1000 and M = 2000 market scenarios, while 
the systematic credit risk factor is discretized with N = 1000, N = 2000 and 
N = 5000 using the method described above. For CVaR calculations, we employ 
the 95% and 99% conhdence levels; This has been dictated by our choice of N. 
As we mentioned earlier, importance sampling methods would be more suitable 
at higher conhdence levels. Note that there would be a dramatic increase in the 
computational cost of simulating the number of exposure scenarios in this case. 

Counterparty exposures (EPEs) are sorted in decreasing order. Let Wn be the largest 
exposure; then the Herfindahl index of the N largest exposures is defined as: 

The effective number of counterparties among the N largest counterparties with respect to total 
portfolio exposures is . 
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Effective Counterparty Number 



Figure 2: Effective number of counterparties for the largest 410 counterparties. 


Effective Counterparty Number 



Figure 3: Effective number of counterparties for the entire portfolio. 


The ranges of individual counterparty exposures are plotted in figure The 95*^ 
and 5*^ percentiles of the exposure distribution are given as a percentage of the 
mean exposure for each counterparty. The volatility of the counterparty exposure 
tends to increase as the the mean exposure of the respective counterparties de¬ 
creases. In other words, counterparties with higher mean exposure tend to be less 
volatile compared to counterparties with lower mean exposure. Given the above 
characteristics, we would expect that wrong-way risk could have an important 
impact on portfolio risk, and that the contribution of idiosyncratic risk will also 
be signihcant. The distribution of the total portfolio exposures from the exposure 
simulation is given in Figure The histogram shows that the portfolio exposure 
distribution is both leptokurtic and highly skewed. It is important to employ 
such highly skewed with very fat tail exposure distribution to ensure the proper 
conservatism of our method in practice. 
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Percentiles of the distribution of exposures for individual counterparties 



120 160 200 240 280 320 360 400 

Counterparty 


Figure 4: 5% and 95% percentiles of the exposure distributions of individual counter¬ 
parties, expressed as a percentage of counterparty mean exposure (Here counterparties 
are sorted in order of decreasing mean exposure). 



Exposure (% base case mean) 

Figure 5: Histogram of total portfolio exposures from the exposure simulation. 


4.2 Numerical Results 


To assess the severity of the worst-case joint distribution, and to determine the de¬ 
gree of conservativeness in earlier methods, we compare risk measures calculated 
using the worst-case joint distribution to those computed based on the stress¬ 


testing algorithm presented in Garcia-Cespedes et ah 2010 and Rosen and Saun¬ 


ders 2010 . In this method, exposure scenarios are sorted in an economically 


meaningful way, and then a two-dimensional copula is applied to simulate the 
joint distribution of exposures (from the discrete distribution dehned by the ex¬ 
posure scenarios) and the systematic credit factor. The algorithm is efficient, and 
preserves the (simulated) joint distribution of the exposures. Here, we apply a 
Gaussian copula, and sort exposure scenarios by the value of total portfolio expo¬ 
sure (this intuitive sorting method has proved to be conservative in many tests. 
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Case I 

MN = 0(10®) 

M = 1000 market scenarios 

N = 1000 credit scenarios 

a 

min^GVaRgys / CVaRwcc) 

max(CVaRsys / GVaRwcc) 

0.95 

52.2% 

95.1% 

0.99 

51.8% 

95.9% 

Case II 

MN = 0(10®) 

M = 2000 market scenarios 

N = 2000 credit scenarios 

a 

min^GVaRgys/CVaRwcc) 

max^CVaRgys j GVaRwcc) 

0.95 

50.3% 

96.9% 

0.99 

49.6% 

96.6% 

Case III 

MN = 0(10'^) 

M = 2000 market scenarios 

= 5000 credit scenarios 

a 

min(CVaRsys/CVaR„cc) 

max(CVaRsyg/CVaRwcc) 

0.95 

44.8% 

96.4% 

0.99 

44.1% 

97.1% 


Table 1: Minimum and maximum of ratio of systematic CVaR using the Gaussian copula algorithm to systematic 
CVaR using the worst-case joint distribution for the largest 220 counterparties at 95% and 99% confidence level. 


Case I 

MN = 0(10®) 

M = 1000 market scenarios 

N = 1000 credit scenarios 

OL 

min(C VaRtot /CVaR„cc) 

max(CVaRtot / CVaRwcc) 

0.95 

52.2% 

96.8% 

0.99 

51.6% 

96.5% 

Case II 

MN = 0(10®) 

M = 2000 market scenarios 

N = 2000 credit scenarios 

OL 

min(C VaRtot /CVaR„cc) 

max(CVaRtot / CVaRwcc) 

0.95 

49.6% 

97.2% 

0.99 

48.9% 

97.4% 

Case III 

MN = 0(10'^) 

M = 2000 market scenarios 

N = 5000 credit scenarios 

OL 

min(C VaRtot /CVaRwcc) 

max(CVaRtot/CVaRwcc) 

0.95 

45.8% 

98.1% 

0.99 

44.7% 

98.2% 


Table 2: Minimum and maximum ratio of CVaR for total losses using the Gaussian copula algorithm to CVaR for 
total losses using the worst-case joint distribution for the largest 410 counterparties at 95% and 99% confidence 
level. 


conducted in Rosen and Saunders 2010| ). For each level of correlation in the 
Gaussian copula, we calculate the ratio of risk (as measured by 95% and 99% 
CVaR) estimated using the sorting method to risk, which in turn, is estimated 
using the worst-case loss distribution. 


We present the results of three discretizations of the worst-case joint distribution. 
Case I employs M = 1,000 market scenarios and N = 1,000 credit scenarios; 
Case II doubles the market and credit scenarios. Lastly in Case III we use 2,000 
market scenarios and 5,000 credit scenarios. Note that Case I and Case II yield a 
discretized worst-case distribution of (9(10®) while Case Ill’s output is of (9(10^). 


Figure shows the graph of the ratio of systematic CVaR using the Gaussian 
copula algorithm described in Rosen and Saunders 2010 to systematic CVaR us¬ 
ing the worst-case joint distribution method described in section 3 for the largest 
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Case I Case II 



correlation 


Case III 

Figure 6: Ratio of systematic CVaR using the Gaussian copula algorithm to systematic 
CVaR using the worst-case joint distribution for the largest 220 counterparties at 95% 
confidence level. 
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Figure 7; Ratio of CVaR for total losses using the Gaussian copula algorithm to CVaR 
for total losses using the worst-case joint distribution for the largest 410 counterparties 
at 99% confidence level. 
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220 counterparties at 95% confidence level. The ratio of the CVaR of the sys¬ 
tematic portfolio loss to the CVaR calculated using worst-case joint distribution 
across various levels of market-credit correlation and discretization scenarios indi¬ 
cate that the distribution simulated using the worst-case joint distribution has a 
higher CVaR compared to previous simulation methods for the largest 220 coun¬ 
terparties by 4.9%, 3.1% and 3.6% at a = 0.95 respectively when the systematic 
risk factor and market risk factor are fully correlated. The difference is larger for 
lower levels of market-credit correlation in the stress testing algorithm. The sort¬ 
ing methods do indeed produce relatively conservative numbers (at high levels of 
market-credit correlation) for this portfolio. In addition to the results presented 
in hgure table [T] shows the minimum and maximum of systematic CVaR to 
CVaR calculated from worst-case joint distribution at 99% conhdence level. Note 
that the results are consistent with what we observed at a = 95%. 

Similar results for calculating CVaR ratios using the total portfolio loss for the 
largest 410 counterparties which constitute more than 99.6% of total portfolio 
exposure are shown in hgure The graphs are based on a higher conhdence 
level, a = 99%, compared to hgure Table shows comparable results to those 
presented in table [T] when we use total portfolio loss instead of systematic loss. 


5 Conclusion and Future Work 


In this paper, we studied the problem of hnding the worst-case joint distribution 
of a set of risk factors given prescribed multivariate marginals and a nonlinear loss 
function. We showed that when the risk measure is CVaR, and the distributions 
are discretized, the problem can be solved conveniently using linear programming. 
The method has applications to any situation where marginals are provided, and 
bounds need to be determined on total portfolio risk. This arises in many hnan- 
cial contexts, including pricing and risk management of exotic options, analysis of 
structured hnance instruments, and aggregation of portfolio risk across risk types. 
Applications to counterparty credit risk were emphasized in this paper, and they 
include assessing wrong-way risk in the credit valuation adjustment, and coun¬ 
terparty credit risk measurement. A detailed application of the algorithm for 
counterparty risk measurement on a real portfolio was subsequently presented 
and discussed. 

The method presented in this paper will be of interest to regulators, who are 
interested in determining how conservative dependence structures estimated (or 
assumed) by risk managers in industry actually are. It will also be of interest 
to risk managers, who can employ it to stress test their assumptions regarding 
dependence in risk measurement calculations. 
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